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^ . Abstract 

o . 

We study the spectrum of H(m) = 75W(— m) with W(m) being the Wilson-Dirac 
i operator on the lattice with bare mass equal to m. The background gauge fields are 

G> \ generated using the SU(3) Wilson action at /3 = 5.7 on an 8^ x 16 lattice. We find 

evidence that the spectrum of H(m) is gapless for 1.02 < m < 2.0, implying that 
the physical quark is massless in this whole region. 
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1 Introduction 



The conventional wisdom obtained from lattice simulations of QCD using 
Wilson fermions at various values of gauge couplings (3 and bare quark mass 
m is that there is a critical line mc{(3) where the physical quark is massless. 
This critical line is not universal, since different definitions for what one means 
by massless quark give different results for mc(/3) at finite (3 because the lattice 
theory has artifacts that only vanish in the continuum limit. One method is 
to define mc{(3) as the line where the pion mass is zero. Such a definition can 
only be implemented by extrapolating from values of m where the pion mass is 
substantially away from zero. In dynamical simulations, critical slowing down 
makes it hard to generate enough independent configurations to obtain a good 
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estimate of the pion mass when it is small. In quenched simulations, there is no 
problem of critical slowing down in the generation of gauge field configurations, 
but it becomes difficult to compute quark propagators below a certain pion 
mass due to the presence of small eigenvalues in the spectrum of the Wilson- 
Dirac operator. Thus, as already mentioned, one relies on an extrapolation to 
estimate mdlS) using data for pion masses greater than O.brrip 

The chiral symmetry of QCD with massless quarks is expected to be sponta- 
neously broken. In the continuum, the spectral density p(A) of 75P is symmet- 
ric, and spontaneous breakdown of chiral symmetry is due to a non-zero value 
for p(0) 0. Furthermore, gauge fields with nontrivial topology will result in 
exact zero eigenvalues of p with a specific chirality 0]. To study the spon- 
taneous breakdown of chiral symmetry on the lattice using Wilson fermions, 
and also to compute the topological charge of the gauge fields using Wilson 
fermions as a probe, we consider H(m) = '~f^W{—m), where W(m) is the stan- 
dard Wilson-Dirac operator with bare mass equal to m. The operator H(m) 
can be viewed as the Hamiltonian describing the evolution of a four dimen- 
sional fermion in a fifth direction. The evolution is in a gauge field background 
that is independent of the fifth direction, and this Hamiltonian plays a cen- 
tral role in the definition of the chiral determinant [Q. Level crossings in the 
spectrum of H(m) as a function of m correspond to the presence of topo- 
logical objects in the background gauge field configuration with the net level 
crossing being the topological charge of the gauge field [||. Level crossings 
at finite values of f3 will not occur at a single value of m, but will instead 
occur over a region in m [^. This is expected to be the case in both pure 
gauge theories and QCD where the fermion dynamics is included. As such, 
the spectrum of H(m) will not have a gap in a region mi(/?) < m < m2{(3). 
In this region, the spectral density in the thermodynamic limit is expected to 
have a finite density of eigenvalues at zero leading to a spontaneous break- 
down of chiral symmetry |0. It is natural to say that the physical quark mass 
is zero in the whole region mi{(3) < m < m2{(3). The existence of a region 
as opposed to a single point is due to finite lattice spacing, and we expect 
both mi{P) and m2(/3) to go to zero as /? goes to infinity. Therefore, one does 
not have a critical line mc{(3) where the physical quark is massless. Instead, 
one has a region in the {P,m) plane bounded by the lines mi(/?) and m2{[3) 
where the physical quark is massless. The critical line mc{(3) obtained from 
the extrapolation of pion masses measured in the region m < mi{f3) will most 
likely lie in the region between mi(/3) and m2(/3), but it need not coincide 
with mi(/3). If mc{(3) is the only point where the physical quark mass is zero, 
then we have a region of bare quark masses between mi{(3) and mc{l3) where 
the physical quark mass is positive but the spectrum does not have a gap. 
Because of low lying eigenvalues in the region mi{f3) < m < mc{(3), it is 
difficult to measure the propagators of mesons. Low lying eigenvalues make 
it numerically difficult to compute the propagator, and even if one succeeds 
in computing the inverse, violent fluctuations in the propagator over the en- 
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semble of gauge fields make it prohibitive to get a good statistical estimate 
of the propagator. Therefore, in the region mi(/3) < m < mc{l3), one cannot 
extract the pion mass by conventional means and show that it is non-zero. 
A method has been recently invented to remove the low lying eigenvalues in 
the region mi(/5) < m < rndP) under the assumption that the low lying 
eigenvalues in the above region are lattice artifacts. This can be thought of as 
a non-local modification of the Wilson-Dirac operator; however, it might en- 
able one to compute the pion mass in the region mi(/3) < m < mc{/3). Gauge 
field configurations with different physical properties are responsible for the 
low lying eigenvalues at different values of m and all these gauge field configu- 
rations in the ensemble are equally important and collectively give rise to the 
existence of the chiral condensate. Therefore, it is not appropriate to remove 
the low lying eigenvalues in the region of mi{(3) < m < mc{(3). Instead, we 
take the point of view that the physical quark is massless in the whole region 
mi(/3) < m < m2{P), which shrinks to the point m = as /3 goes to infinity. 

Our aim in this paper is to present some evidence in support of the picture 
outlined in the previous paragraph. To this end, we focus on pure gauge field 
configurations generated using the SU(3) Wilson action at /5 = 5.7 on an 
8^ X 16 lattice. Using an ensemble of 50 configurations, we will show that 
H(m) is gapless in the mass range between mi = 1.02 and 1712 = 2.0, indi- 
cating that the physical quark is massless in this whole region. We do not 
consider m > 2 since this is a region of m which is clearly unphysical in the 
continuum limit. Using the method proposed in Ref 0, we will show that 
chiral symmetry is spontaneously broken in this whole region of m. For the 
quenched approximation on a 16^ X 32 lattice at /3 = 5.7, the data of Ref. M 
at m < 1 yield a an extrapolated value of rric = 1.047 where the pion mass 
becomes zero. This is inside the gapless region. We will show that the pion 
propagator for a given gauge field background can be estimated quite well 
in this region by using the spectral representation for H(m) and restricting 
the spectral sum to the low lying eigenvalues. This will enable us to compute 
the pion propagator more efficiently than with conventional means since we 
can sum over all possible source points on the lattice and effectively restore 
translational invariance. The resulting propagator, when averaged over an en- 
semble of configurations, will still have large fluctuations due to the variation 
in the magnitude of the low lying eigenvalues. But close to mi = 1.02 the 
fluctuations will not be large, and we will see that the pion propagator in the 
region from mi to m^ is not described by an exponential decay but rather by a 
power law, providing further evidence that the physical quark described by the 
Wilson-Dirac fermion is massless in this region. A little more physical insight 



In fact, m > 1 does not yield a unitary theory in the continuum limit, but we 
are, however, away from the continuum Umit here. At m = 2 the free fermion has 
four poles in the propagator due to "doublers" , and since we are interested in the 
gap in the spectrum of the H(m) we avoid the region of m > 2. 
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on the properties of the gauge fields responsible for the low lying eigenvalues 
will be gained by looking at the eigenvectors of H(m) corresponding to these 
low lying eigenvalues. In particular, the eigenvectors corresponding to zero 
eigenvalues of H(m) have a natural association with localized objects on the 
lattice, and the eigenvectors themselves should be localized. In the sample of 
50 configurations, we have nine such zero eigenvalues of H(m) for m < 1.047. 
We will show that the corresponding eigenvectors are all localized. The zero 
spatial momentum component of these eigenvectors has a large extent in time, 
occupying anywhere from seven to ten lattice spacings. A detailed study of 
the level crossing over the whole range of m from m = 1 to m = 2 yields in- 
formation about the topological content of the gauge field, including the size 
distribution of the topological objects on the lattice. The analysis of the size 
distribution shows that larger objects cross at smaller values of m, indicating 
that the crossings at m < should indeed correspond to larger objects than 
the ones associated with crossings at m > m^. 

This paper is organized as follows. In the next section we review the connection 
between the level crossings of II(m) and topology. We also review the method 
to measure the chiral condensate using a term that breaks parity This 
method measures the chiral condensate arsing from spontaneous breaking only 
and gives a zero result unless the spectrum is gap less. In section 3, we measure 
the chiral condensate using the spectral representation of H(m) and show that 
there is a non-zero chiral condensate in the whole region from mi = 1.02 to 
7712 = 2.0. We also obtain the spectral density as a function of m and show 
that the gap closes at mi = 1.02, and we find that it remains closed for all 
higher values of m up to the largest value we studied, namely m2 = 2. This is 
consistent with a non- vanishing chiral condensate in this whole region of m. 
In section 4, we focus on the region from mi = 1.02 to rric = 1.047. In this 
region we compare the shape of the eigenvectors very close to the crossings 
points with the pion propagators at the same value of m. We find excellent 
agreement, indicating that the pion propagator is dominated by that particular 
eigenvector. By employing the spectral representation of H(m) and taking a 
restricted spectral sum, we show that one can obtain the pion propagator 
quite well over a region around the crossing point. We then use the restricted 
spectral sum to compute the average pion propagator over the entire ensemble 
and show that the resulting pion propagator is not described by an exponential 
decay. Our conclusions are presented in section 5. 
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2 Properties of the spectral flow of H(m) 



We start this section by writing down the central quantity in the paper, namely 
the Wilson-Dirac Hermitian Hamiltonian H(m) on the lattice in a fixed gauge 
field background. Upon suppressing the dependence on the SU(3) background 
gauge field U, H(m) in the chiral basis is 

HM=(°Ct'" -B + ™)' (1) 
where 



B,o.,Ak, k') = l6o.,pY.['^6i,6kk' - U;^{k)6k',k+f. - {\JlTik')b^^y^(:\ (3) 

^^ = (? J)' ^^=(° 7)' ^^=(J -l)' ^^=(o z) ' 

and where and k' label the lattice sites, a and /5 are two component spinor 
indices, and i and j are three component color indices. 

The fermionic action in terms of H(m) is 

= iP'R{m)ilj, (5) 

with ip' related to the standard ip hj ijj = 4^'^^. The fields ip' and ip are 
independent Grassmann variables. The choice for the sign of m differs from 
the standard one. In the free theory {U = 1), the action with — m describes a 
free fermion with mass equal to m. 

The operator H(m) can also be viewed as the single particle Hamiltonian 
describing the propagation of a fermionic field in the fifth direction under 
a background gauge field that is static in that direction. The many-body 
Hamiltonian describing the fermionic evolution in the fifth direction is 

n{m) = a^R{m)a, (6) 

where a and are canonical fermion annihilation and creation operators for 
Dirac fermions. The chiral determinant on the lattice, which is the result of 
integrating out a single Weyl fermion in a fixed gauge background, is given by 
the overlap |Q 

detC = (0- |0+) , (7) 
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where |0±) are the ground states of T-C{±m) for some fixed value of m in the 
range (0, 2). All choices of m in this range give rise to the same chiral theory 
in the continuum limit but will have different lattice spacing effects away from 
the continuum. 

Equation (|^) along with provides an immediate connection between level 
crossings in the spectral flow of H(m) as a function of m and the global 
topology of the background gauge field On a finite lattice containing V 
sites, H(m) is a 2K x 2K matrix with K = 6V. The ground states |0±) are 
obtained by filling all the negative energy states of H(±m). For an arbitrary 
gauge field background, one can rigorously show that H(— m) always has K 
negative energy states for all values of m > 0. Consider a background gauge 
field for which H(m) has {K — Q) negative eigenvalues with Q > 0. Then det C 
as defined in is zero and {0 — \a\ ■ ■ ■ aQ\0+) is non-zero, implying that Q 
fermions of a fixed chirality are created in this gauge background. Such a gauge 
field configuration on the lattice will then be said to carry a topological charge 
of Q in analogy with the continuum index theorem [Q. Examples of such gauge 
field configurations can be constructed on the lattice 0]. Since H(— m) has K 
negative energy states, while H(m) has {K — Q) negative energy states, the 
spectral flow of H(m) as a function of m should have a net level crossing of 
Q (from below zero to above zero) as one goes from —m to m. Therefore 
H(m) should have at least Q zero eigenvalues, and these will generically occur 
at Q different values of m. Typically, n+ levels will cross from below zero to 
above zero, and n_ levels will cross from above zero to below zero, with the 
topological charge being given hy Q = n^—n_. Then H(m) will have (n++n_) 
zero eigenvalues at (n_(. + n_) different values of m. 

It is instructive to write down the zero eigenvalue equation of H(m) using (|l]): 

(B - m)u + Ct; = 0; C^u - {B - m)v = u^u + v'^v = 1. (8) 

Trivial manipulation of the above equation results in the following necessary 
condition for a zero eigenvalue [Q 

m^Bm + f ^Bf = m. (9) 

Since B, defined in (^, is positive for all values of the background gauge field, 
it follows that zero eigenvalues occur only for m > @] • As one takes a lattice 
gauge field and interpolates to the continuum, all solutions to (§) that occurred 
at some value of m > will now occur at m = 0, since the norm of B on any 
state goes to zero in the continuum. Level crossings on the lattice will occur 
at different values of m because the eigenvectors associated with the level 
crossings have a typical size inside which they are localized. Eigenvectors with 
larger sizes will vary slower across the lattice than eigenvectors with smaller 
sizes, implying that m as given in (P) will be bigger for smaller size objects. 
Since we expect topological objects of varying sizes in an ensemble of gauge 
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fields generated at some finite gauge coupling on the lattice, we expect level 
crossings to occur in a region of m. 

Gauge fields generated at a finite gauge coupling on the lattice have a finite 
correlation length in lattice units. One can define the spectral distribution 
of H(m) in the thermodynamic limit, in which the number of lattice sites 
becomes large, as the one obtained by considering a large ensemble of gauge 
field configurations generated on a lattice with all linear sizes greater than the 
correlation length. The spectrum of H(m) corresponding to each gauge field 
configuration in the ensemble will be discrete, but the spectral distribution for 
an infinite ensemble is expected to be continuous. We denote the continuous 
spectral distribution by p{X;m). This will be an even function of A since for 
every eigenvalue A of H(m) in a fixed gauge field configuration there is an 
eigenvalue —A for the Hamiltonian associated with the parity transformed 
gauge field on the latticei. If p(0; m) is non-zero, then H(m) does not have 
a gap at that value of m. Low lying eigenvalues of H(m) and level crossings 
in H(m) can give rise to a non-zero value for p(0;m). One can establish a 
non-zero chiral condensate if p(0; m) is non-zero by adding a term ihip'tlj to 
the fermionic action in (]^) In the presence of this term, 

in a fixed gauge field background U, where 2K is the total number of degrees 
of freedom. The right hand side is obtained by taking the average over U and 
its parity transformed partner. To obtain the chiral condensate, one has to 
compute the ensemble average in the limit of infinite volume at finite values 
of h and extrapolate to the limit of vanishing h. In this limit, we obtain a 
connection between the chiral condensate and the non-zero value of p(0; m) 0, 
namely, 

i{iP'^) =7rp{0;m) . (11) 

Note that a non-zero value for comes solely from p(0;m). This is the 

advantage of measuring the chiral condensate using Existence of the 

chiral condensate can be obtained by an explicit computation of (|10D over an 
ensemble of configurations function of h. 



3 Evidence for a chiral condensate over a region of m 

We generated 50 independent SU(3) gauge field configurations on an 8^ x 16 
lattice at /3 = 5.7 using the standard single plaquette Wilson action. We 

^ Lemma 4.4 in Ref [Q] provides a proof of this statement. 
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Fig. 1. Spectral distribution of H(m-) at m = 1 obtained from the ten eigenvalues closest 
to zero on 50 gauge field configurations at /3 = 5.7. The bin width is 0.005. 

computed the ten eigenvalues of H(m) closest to zero as a function of m. For 
this, we used the Ritz functional method developed in Ref. |^ and improved 
for the case of multiple eigenvectors in Ref. [|T0[. The Ritz functional method 
finds the lowest eigenvalues of H^(m), and these are in fact the eigenvalues of 
H(m) closest to zero. The ten eigenvalues of H^(m) were obtained to a relative 
precision of 10~^ or to an absolute precision of 10~^ depending on which 
condition is satisfied first in the iterative procedure. Only ten eigenvectors were 
kept in the iterative procedure and were reorthogonalized every ten iterations. 
Eigenvalues and eigenvectors were obtained for masses between m = 1 to 
m = 2 with increments of Am = 0.02. In first order perturbation theory 

= -0t,M750nM , (12) 

with 0„(m) being the eigenvector corresponding to A„(m). The mass increment 
was sufficiently small to allow for an accurate interpolation of eigenvalues 
between consecutive values of m. The ten eigenvalues closest to zero varied in 
the range from —0.1 to 0.1 as a function of m over the entire ensemble. Figure 
shows the spectral distribution p(A; m) at m = 1. This distribution has a gap 
and the physical quark is massive. Figure ^ shows the spectral distribution 
p(A; m) at m = 1.03. Here the spectrum does not have a gap and the physical 
quark is massless. 

Since p(0; 1.03) is non-zero, we expect a chiral condensate at m = 1.03. To 
investigate the chiral condensate as a function of m, we use (0) to measure 
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Fig. 2. Spectral distribution of H(m) at m = 1.03 obtained from the ten eigenvalues 
closest to zero on 50 gauge field configurations at (5 = 5.7. The bin width is 0.005. 
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Fig. 3. The chiral condensate i{il)'il)){m) from m = 1.0 to m = 1.1 with h = 0.001. The 

middle curve is the estimate of the chiral condensate, and the curves on cither side show 
the error of this estimate. The solid line at the bottom is the result for free fermions. 



9 



0.004 

? 

A 

1^ 0.002 



0.000 

1.0 1.2 1.4 1.6 1.8 2.0 
m 

Fig. 4. The chiral condensate ijj){m) from m = 1.0 to m = 2.0 with h = 0.001. The 
middle curve is the estimate of the chiral condensate, and the curves on either side show 
the error of this estimate. The solid line at the bottom is the result for free fermions. 

i{ip'ip)u per configuration and then compute an ensemble average to obtain 
For small values of h in (|T0|) , only the low lying eigenvalues will con- 
tribute. Setting h = 0.001 in ( ]TI1| ) and summing over the ten lowest eigenvalues 
of H(m), we obtain a value for z (■?/''?/') (m) as a function of m. We also compute 
i{ijj'ip){m) for free fermions ai h = 0.001 by summing over the whole spectrum. 
We plot the result in the range 1 < m < 1.1 in Fig. ^. The solid line at the 
bottom of the figure is the result for free fermions. When m > 1.02, the result 
from the gauge theory is significantly larger than from free fermions, and is 
consistent with a non-zero chiral condensate in this mass region. Figure. § illus- 
trates the behavior of the condensate i{ip'ip){m) for the region 1.0 < m < 2.0. 
This result is consistent with a non-zero value for p(0;m) for the entire range 
of masses between mi = 1.02 and m2 = 2.0. Note that the critical quark mass 
rric = 1.047, obtained by extrapolating the pion mass to zero, is inside the 
region where the spectrum is gapless. To understand the dependence of the 
chiral condensate on the parameter h, in Fig. ^ we plot versus /i at a 

typical mass m = 1.5. The picture is basically the same at other values of m 
greated than mi. Clearly there is a wide region of h where the condensate is 
basically flat, in contrast to the free fermion result which drops exponentially 
in the plot with a logarithmic scale. As expected from (PUj), the condensate 
goes to zero as h goes to zero and as h goes to infinity. The errors at the peak 
region are large due to the large fluctuations in coming from the low lying 
eigenvalues. If the spectral sum is restricted to the six lowest eigenvalues, as 
opposed to the ten lowest, all the above results are essentially unaltered. 
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Fig. 5. Chiral condensate i{ip'Tp){h) at m = 1.5 viewed as a function of h. The x-axis is 
logarithmic to show a wide range in h. The middle curve with a peak is the estimate of 
the chiral condensate, and the curves on either side show the error of this estimate. The 
exponentially varying solid line is the result for free fcrmions. 

4 The mass range mi < m < rric 

In the previous section we explored the consequences of the low lying eigenval- 
ues of H(m). We showed that the spectrum is gapless in a mass range bounded 
by nil = 1-02 and m2 = 2.0, and we provided evidence for a non-vanishing 
chiral condensate in the gapless region. Our interpretation is that the physi- 
cal quark is massless in this whole region of m. If this is indeed the case, we 
should find zero eigenvalues of H(m) scattered throughout this region due to 
the presence of topological objects in the gauge field configurations. Equation 
(H) defines the zero eigenvalue equation for H(m), and is a necessary condi- 
tion for a zero eigenvalue. Based on (P), it was argued in section 2 that larger 
topological objects will result in a zero eigenvalue at smaller values of m. In 
this section we first look at the sizes of the topological objects and their cross- 
ing points, and we present evidence supporting the above statement. Then we 
explore the region mi < m < rric, where rric = 1.047 is the point at which the 
extrapolated pion mass vanishes, to establish the following points: 

- Eigenvectors associated with zero eigenvalues of H(m) in the mi < m < m^ 
region correspond to large localized objects extending roughly over half the 
lattice in the long direction. 
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- The pion propagator for a fixed gauge field background with level crossings, 
at values of m close to a crossing point, is well described by a spectral sum 
restricted to a few low lying eigenvalues of H(m). 

- The ensemble average of the pion propagator obtained by a restricted spec- 
tral sum at mi = 1.02 do not fit an exponential decay. Instead they fit a 
power law. 

To this end we write down the spectral representations of H(m) and H^^(m). 
The eigenvalue equation for H(m) is 

Y^R,,,,,{k,k')<p-^{k') = x^ram , (13) 

k'bj 

and we have dropped the dependence on m in the above equation. Here k and 
k' label the sites on the lattice, a and b are four component spinor indices, 
i and j are three component color indices, and n labels the eigenvalues. The 
eigenvectors are normahzed to yield Z]fcaj0S(^)]*0ai(^) = ^nm- In terms of the 
eigenvalues and eigenvectors, the spectral representations of H and are 

Ha.,6,(A:, A;') = E ^n0:.(fc)[0^,(A:')]* (14) 

n 

Kl,Ak,k') = C.(^)[C,(^')]* • (15) 

The point-like pion field is given by 

vr(A:)=EC(fc)^«(fc) , (16) 

ai 

and using the spectral representation of H^^ in ([TSD , the pion propagator in a 
fixed gauge field background becomes 

{7r{k)^{k'))u = E Y^fnmik)f:mik') , (17) 

where the gauge invariant functions fnmik) are defined by 

fnm{k) = E[C(fc)]*C(fc) • (18) 

ai 

The eigenvalues of H{m) will sweep out a spectral fiow as m varies over its 
range, and typically some levels will cross, thereby producing a set of zero 
modes at discrete values of the mass. Let us examine such a zero mode, 0aj(fc), 
at some crossing point m. In what follows, it is convenient to consider the zero 
momentum part of /oo(^) defined by f{t) = J2k foo{k,t). We can associate a 
size with f{t) using the pion correlator 

^o(t) = E/(0/(t + , (19) 
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Fig. 6. The localization size as a function of the crossing point. The data point at m includes 
all crossings that occurred in the interval m it 0.05. 

which is a periodic function of t with a peak at t = 0. We define the locahzation 
size, p, through the second moment of Co{t): 

Recall that the lattice has 16 sites in the time direction, and in calculating 
the localization size we order them to run from —7 to 8. Figure ^ shows the 
localization size of the zero modes as a function of the location of the crossing 
point. The trend in this figure is consistent with the picture that localized 
objects of a larger size cross at smaller values of m. 

We now focus on crossings in the region between mi and mc, i. e. the mass range 
1.02 < m < 1.047. In our ensemble of 50 configurations, we found seven which 
contained one crossing each within this mass range, and one configuration with 
two crossings between mi and m^. For all nine zero modes, f{t) peaks at some 
value of t in the range < t < 15, and in Fig. |^ we have plotted f{t) for two of 
them. In numerically solving for the eigenvalues of H{m), we generically find 
near zero-modes since we sweep over discrete values of the mass parameter. 
The solid line in Fig. ^ corresponds to a mode whose eigenvalue was closest to 
zero along the computationally obtained discrete spectral flow. It occurred at 
a mass of m = 1.04, and using (|12]) we find that the associated zero crossing 
happens at mo = 1.0375, which is very close indeed. As one can see, the 
mode peaks somewhere between t = 8 and t = 9 and has a width of about 
seven or eight lattice spacings. The dashed line is the eigenvector for another 
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Fig. 7. The function fit) for two near zero-modes within the mass range mi < m < rric- 
The sohd Hne corresponds to a mode at m = 1.04 for a certain configuration within the 
ensemble, while the dashed line corresponds to a mode at m = 1.02 for a different gauge 
field configuration. 

configuration in the ensemble whose eigenvalue along the discrete spectral 
flow is again closest to zero. It occurs at m = 1.02, and ([I^ ) implies that 
the associated zero crossing is at mo = 1.025. This mode has a double peak 
structure, and it stretches over ten lattice spacings. Of the nine zero modes 
within the range mi < m < rric, seven of them had a single peak in f{t) while 
two exhibited double peaks, and all them stretched over seven to ten lattice 
spacings. 

We also looked at foo{k) over the whole lattice and found that the modes were 
localized in all four directions but were not spherical. These modes are large 
objects having widths of seven to ten lattice spacings in the long direction and 
widths of six or so in the short directions. The flnite size of the lattice plays a 
significant role in the shape of the mode. Furthermore, the associated lattice 
gauge field backgrounds of these eigenvectors are far from the continuum limit, 
and as such are not smooth configurations. This also affects the shape of these 
modes. 

Near values of m where H(m) has a zero eigenvalue, we expect the spectral 
decomposition of H~^(m), as given by (|l^), to be dominated by the low lying 
modes. We test this on the configuration corresponding to the solid line in 
Fig. ^ at two mass points m = 1.029 and m = 1.036, the latter being close to 
the crossing point mo = 1.0375. We shall restrict the spectral sum in (0) to 
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Fig. 8. Comparison between tlie fnll pion propagator (dashed) and the pion propagator 
obtained from the restricted spectral sum (solid), both of which are in the fixed gauge field 
background corresponding to the solid line of Fig. 7. The mass is m = 1.029, which is near 
the crossing value mo = 1.0375. 
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Fig. 9. Comparison between the full pion propagator (dashed) and the pion propagator 
obtained from the restricted spectral sum (solid) using the same gauge field background 
as in Fig. 8. This time the mass is m = 1.036, which is even closer to the crossing value 
mo = 1.0375, and the two curves lie almost on top of one another. 
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run only over the ten eigenvalues closest to zero, and we compute the quantity 
C{t) = Z]fc(7r(0)7r(fc, using both the restricted spectral sum and the full 
propagator obtained by inverting H(m) with standard techniques. Fig. ^ shows 
C{t) at m = 1.029 and Fig. ^ illustrates a comparison at m = 1.036. At 
m = 1.029 both methods agree for 3 < t < 12, while at m = 1.036 they agree 
at all values of t, which is an expected result since the latter mass is closer to 
the crossing point. Furthermore, it should come as no surprise that the pion 
propagator itself, C{t), looks very much like f{t) for the corresponding zero 
mode, in this case the solid line of Fig. [7|. 

The pion propagators of Figs. ^ and ^ are far away from an expected exponen- 
tial decay starting at t = 0, and instead they have a dominant peak around 
t = 8. Such configurations are the reason why one cannot measure the pion 
mass at values of m in the region between mi = 1.02 and rric = 1.047, since 
they cause large statistical fluctuations. This has been known for quite some 



time and was first observed in Ref. |11|. As previously mentioned, the peak is 



due to a zero of H(m) close to values of m where the pion propagator is being 
measured. Since we have presented evidence that the pion propagator can be 
well approximated by a restricted spectral sum, one can reduce the statistical 
fluctuations by computing the pion propagator as an average over all possi- 
ble translations of the background gauge field. This can easily be done in the 
restricted spectral sum, but it would be prohibitive if one were using conven- 
tional methods to compute the pion propagator, since one would have to do 
a volume worth of inversions of the Wilson-Dirac operator. It should be men- 
tioned that a restricted spectral sum has been used to study the spectroscopy 
in the region where H(m) has a gap []T2[. The average pion propagator over 



all possible translations of a fixed gauge field background is 

aAk) = ^j:{^{k')n{k + k')) , (21) 
^ k' 



and we find it useful to study C'{t) = J2kCav{k,t). In Fig. [TO we plot C'{t) 



at m = 1.029 for the same configuration as in Fig. p. This function has the 
expected behavior for a pion correlator in the sense that it falls off starting at 
t = 0, but it does not fit an exponential decay. Instead, it fits a power law of 
the form obtained using the zero mode solution in the presence of a classical 
instanton. The spherically symmetric zero mode associated with a classical 



instanton take the form [jT3| 



i^^xmx) ~ , (22) 

which gives a contribution to the pion correlator of 

^^'^^^ " / (2^ (^^TTF ((^+^F+7F ■ ^^^^ 
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Fig. 10. The solid line is the pion propagator of Fig. 8 averaged over all possible translations 

of the gauge field. The dotted line that falls on top of the solid line is a fit to a power law, 
and the dashed line shows an exponential fall off in comparison. 




Fig. 11. The data points with error bars show the pion correlator C'(t) averaged over the 
ensemble of 50 configurations. The dotted line is a power law decay and the dashed line 
shows an exponential fall oS^ in comparison. 
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The above integral can be evaluated and the resulting functional form can be 
fitted to C'{t). The fit corresponding to p = 5.1375 is shown as the dotted line 
in Fig. [10], which falls right on top of C'{t) represented by the solid line. In 
contrast to this fit, the dashed line shows the behavior of C(e~"* + e""'-^^"*-') 
with a = 0.3. Note that C'{t) is below the exponential decay near t = im- 
plying that the decay is slower at t = when compared to the decay at t = 8. 
This is atypical behavior for the correlator since excited states would result in 
a faster decay at t = 0. This is another indication that the behavior is incon- 
sistent with an exponential decay. The power law behavior is consistent with 
a non- vanishing chiral condensate in the region 1.02 < m < 1.047, as found 
in the previous section. Furthermore, Fig. |^ shows that the zero eigenvalues 
in the above region of m are caused by large topological objects extending 
over several lattice spacings. Therefore, zero eigenvalues are caused by phys- 
ical objects in the background gauge field, which results in a non-vanishing 
chiral condensate and a power law behavior in the pion mass. The proposal 
of Ref. 10] shifts the position of these zeros to some value of m greater than 
rric = 1.047, thereby removing all low lying eigenvalues in the region m < rric. 
This should also make the chiral condensate vanish for m < rric and remove 
the power law contribution to the pion propagator in this region. This indi- 
cates that the proposal in Ref. is a non-local change to the Wilson-Dirac 
operator motivated by the need for a single value of m where the pion mass 
vanishes. The proposal of Ref. 0] does nothing to the zero crossings that occur 
at m > rric, and therefore one will still have a region where H(m) is gapless. In 
this region there will be a non-vanishing chiral condensate and the pion will 
be massless. In effect, the modification of the eigenvalues proposed in Ref. |0] 
slightly changes the region where the spectrum of H(m) is gapless and we see 
no compelling reason to make this non-local modification to the Wilson-Dirac 
operator. 



We now proceed to compute C'{t) as an average over the whole ensemble at 
rrii = 1.02. At this mass value there is evidence for a non- zero density of 
eigenvalues at zero, as presented is section 3. Figure |Tl] shows the result of 
averaging over the whole ensemble. The errors, although large, still enable us 
reasonably estimate of the pion correlator at this value of m. The dotted line 
shows the same power law function as the one used in Fig. |10|with p = 5.675, 
whereas the dashed line shows an exponential decay with a = 0.25. Again we 
see the same behavior as in Fig. |10|, namely the data is below the exponential 
decay at t = 0. The pion propagator at different times is correlated and 
it is interesting to note that the power law fit goes through all the central 
values. The propagator fits the power law better and is not consistent with 
an exponential decay. We have therefore shown that it is not appropriate to 
alter the low lying eigenvalues in the region of m between mi = 1.02 and 
rric = 1.047. We have also established that (|21| ) is a powerful method for 
estimating the pion propagator in this region using the restricted spectral 
sum to compute the right hand side of this equation. 
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5 Conclusions 



We have studied tlie spectrum of H(m) = 75W(— m), witli W(m) being the 
Wilson-Dirac operator on the lattice with a bare mass m. We generated 50 
pure SU{3) gauge configurations using the single plaquette Wilson action at 
P = 5.7 on an 8^ x 16 lattice. We found that the spectrum is consistent 
with it being gapless in the region between mi = 1.02 and 1712 = 2.0. Using 
a parity breaking operator as proposed in Ref. we estimated the chiral 
condensate and found it to be non-zero in this whole region. Therefore, it 
is natural to say that the physical quark is massless in this region. This is 
in contrast to declaring a single value, rric = 1.047, where the pion mass 
obtained by extrapolation from the region m < 1 vanishes f^, as the unique 
point where the physical quark is massless. We computed the pion correlator 
using a restricted spectral sum and by averaging over all possible translations 
of the background gauge field. We found that the pion correlator is consistent 
with a power law decay at mi = 1.02, with the power law coming from the 
low lying eigenvalues of H(m). This further confirms that the physical quark 
is massless in this whole region of m, as opposed to a single point. We also 
showed that the eigenvectors corresponding to zero eigenvalues of H(m) in the 
region between mi = 1.02 and rric = 1.047 are localized and roughly span 
half the lattice in the long direction (about 1.36 fm). As such, the low lying 
eigenvalues in this region are associated with physical objects in the gauge 
field background. 



The region in m where the physical quark is massless shrinks to the single point 
m = in the continuum limit. But the region has a finite width if j3 is finite. 
The width will change if we improve the Wilson-Dirac operator. We studied 
the spectral fiow for a few gauge field configurations using the Sheikholeslami- 
Wohlert improvement scheme. We found that the gap closes at a smaller 
value of m, but preliminary indications are that it remains closed all the way 
up to m = 2, as in the Wilson-Dirac case. For the Wilson-Dirac operator, 
the width of the gapless region is not expected to change as we go to infinite 
lattice volume at fixed (3. But we expect more low lying eigenvalues as we go 
to larger lattice volume since the chiral condensate arises from spontaneous 
symmetry breaking, and the effect should become stronger in larger lattice 
volume. For the same reason, one should see more zero eigenvalues of H(m) 
in a fixed region of m. An investigation of ten configurations on a 16^ x 32 
lattice at /3 = 5.7 showed that there were seven zero eigenvalues in the region 
1.02 < m < 1.047. This is a greater percentage than the nine zero eigenvalues 
found in 50 configurations on the 8^ x 16 lattice at the same coupling. The 
eigenvectors associated with these modes were found to be locahzed and large 
just like on the 8^ x 16 lattice. There will be two competing effects as one goes 
to larger volume. As remarked, one will find more zero modes in the region 
1.02 < m < 1.047. But the effect of a single zero mode in larger volume will 
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not be felt in as wide a region around the crossing point as in a smaller volume. 
This is because the number of modes that contribute to the pion propagator 
increase with volume and each single mode therefor has less weight. 

Since for /3 = 5.7 the gap does not reopen, the natural conclusion is that we are 
far from the continuum limit. We expect the gap to open again before m = 2 
if we go to a larger (3. Tracing out the gapless region in the (/?, m) plane 
is necessary to probe massless quarks as one goes to the continuum limit. 
It is also important for a proper study of massless QCD using the overlap 
formalism P| or the domain wall formalism |T3]. For a proper study one has 



to work at a value of (3 so that the gapless region ends before m = 2. To 
reproduce the proper chiral behavior, the value of the mass that goes into the 
overlap Hamiltonian in (^, and the domain wall mass, should be kept greater 
than the upper bound of the gap m2{[3). The gap also has to open before 
m = 2 so that one can use level crossings for a proper measurement of the 
topological charge and topological susceptibility. 
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